arXiv: math. PR/0000000 



Bandwidth sharing networks 
with priority scaling 

M. Feuillet, M. Jonckheere and B. J. Prabhu 

INRIA Paris- Rocquencourt, 
Domaine de Voluceau, Rocquencourt, B.P. 105 
78153 Ic Chesnay Cedex, France, 
e-mail: mathieu . f euillet{Binr ia. f r 

CONICET, Departamento de Matemdtica, 
Facultad de Ciencias Exactas y Naturales, 
Universidad de Buenos Aires, 
Pabellon 1, Ciudad Universitaria, 
1428 Buenos Aires, Argentina, 
e-mail: mj onckheOdm . uba . ar 

CNRS; LAAS; 
7 avenue du Colonel Roche, 
F-31077 Toulouse, France. 
Universite de Toulouse; 
UPS, INSA, IMP, ISAE; LAAS; 
F-31077 Toulouse, France, 
e-mail: balaOlaas.fr 

Abstract: In multi-class communication networks, traffic surges due to 
one class of users can significantly degrade the performance for other classes. 
During these transient periods, it is thus of crucial importance to imple- 
ment priority mechanisms that conserve the quality of service experienced 
by the affected classes, while ensuring that the temporarily unstable class 
is not entirely neglected. In this paper, we examine the complex interaction 
occurring between several classes of traffic when classes obtain bandwidth 
proportionally to their incoming traffic. We characterize the evolution of 
the network from the moment the initial surge takes place until the sys- 
tem reaches its equilibrium. Using an appropriate scaling, we show that the 
trajectories of the temporarily unstable class can be described by a differen- 
tial equation, while those of the stable classes retain their stochastic nature. 
A stochastic averaging phenomenon occurs and the dynamics of the tem- 
porarily unstable and the stable classes continue to influence one another. 
We further proceed to characterize the obtained differential equations and 
the stability region under this scaling for monotone networks. We illustrate 
these result on several toy examples and we finally build a penalization rule 
using these results for a network integrating streaming and elastic traffic. 
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1. Introduction 

Communication networks are dealing with heterogeneous sources of traffic hav- 
ing different behaviors in terms of volume of data and aggressivity. Ideally, the 
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network should respond to the different demands in the fairest possible way, 
i.e. by avoiding a significant degradation of quality of service of a given class of 
traffic when another class undergoes a major traffic surge. 

The impact of large-scale traffic surges, also known as flash-crowds, on web 
servers and content distribution networks has been the subject of several studies 
(Stavrou, Rubenstein and Sahu, 2004; Kandula et al., 2005; Deshpande et al., 
2007). These mainly focus on designing mechanisms to make the content providers 
resilient to surges of a given type of traffic. However, in addition to overloading 
the content providers, a traffic surge can also negatively impact the performance 
of concurrent flows in the network. The temporarily unstable class can poten- 
tially starve the other classes from network capacity thereby subjecting them 
to unreasonable delays and packet losses. In such circumstances, in addition to 
protection mechanisms in web servers, it is crucial to implement bandwidth- 
sharing mechanisms inside the network that would protect the stable classes 
from the adversarial effects of the surge. It seems natural that such mechanisms 
should penalize the temporarily unstable class more when the level of congestion 
it creates is larger, without actually throttling it. (Thus, the more significant 
the surge is, the smaller the bandwidth each flow in this class gets.) The eon- 
sequences of traffic surges on the performance of the different classes in the 
presence of such bandwidth sharing mechanisms have not been explored much. 

In this paper, we take a global view at the effects of a traffic surge in a multi- 
class communication network: our aim is to present an analytic treatment of 
the complex interaction that takes place between the temporarily unstable class 
and the stable class during a traffic surge when the temporarily unstable class 
is penalized proportionally to its level of congestion. 

Towards this end, we consider stochastic networks describing the evolution of 
the number of flows in a network where different classes of traffic compete for the 
bandwidth. Bandwidth-sharing network models (Massoulie and Roberts, 2002; 
Bonald and Prouticre, 2003; GromoU and Williams, 2009) have become quite 
a standard modeling tool over the past decade. In particular, they have been 
used extensively to represent the flow level dynamics of data traffic in wireline or 
wireless networks (Bonald et al., 2006), as well as for the integration of voice and 
data traffic (Bonald and Proutiere, 2004), hence generalizing more traditional 
voice traffic models, e.g. (Kelly, 1979). 

To obtain structural results, we introduce a scaling when possibly only a 
subset of classes have initial conditions converging to infinity. Those classes 
shall be the ones undergoing a traffic surge (surging classes). We consider a 
situation where the allocation of bandwidth shall be meanwhile weighted such 
that the other classes of traffic (stable classes) are not led to starvation, i.e., 
the priority weight is very small compared to the offered traffic. Accelerating 
time together with re-scaling the state of the surging classes allow to "zoom 
out" the process, just as for usual fluid limits and obtain a bird's-eye view of 
the large scale dynamics for these classes. In order to obtain a classical fluid 
limit for Jackson networks (Robert, 2003) or for more complex bandwidth- 
sharing networks (Gromoll and Williams, 2009), all the classes are jointly scaled 
in time and in space. This yields a set of differential equations that govern 
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the dynamics of all the classes. Under additional assumptions on the drift S 
of the considered Markov process, the differential equation is simply of the 
form x{t) = 5{x{t)) (see the considerable amount of work on fluid limits and 
ODE methods both for Markov processes and for communications networks 
(Dai, 1995; Darling and Norris, 2008; GromoU and Wilhams, 2009; Meyn, 2008; 
Robert, 2003)). 

In our case, the situation differs as the transitions of surging classes arc also 
scaled to model that the priority weight of surging classes is inversely propor- 
tional to the level congestion. This has far-reaching consequences for the struc- 
ture of the limiting process. Under this scaling, we will show that the dynamics 
of the surging classes can be described by a set of deterministic differential equa- 
tions, while the stable classes retain their stochastic nature. Hence, a time-scale 
separation occurs: the surging classes evolve on the much slower time-scale com- 
pared to the stable classes. However even with this separation of time-scales, a 
strong coupling in the dynamics of the surging classes and the stable classes re- 
mains. The dynamics of the surging classes are influenced by the stable classes 
through their conditional distribution which in turn depends on the level of 
congestion of surging-class flows being fixed to their present macroscopic value. 
Hence, for the surging classes the differential equations obtained are of the form 
x{t) = S*{x{t)), where S* is an average of the first coordinates drift according to 
the conditional distribution of the other classes, given the state of the surging 
classes. This phenomenon is usually known in the probability literature as aver- 
aging principle and has been studied by several authors. We follow in particular 
the methodology introduced in the seminal paper (Kurtz, 1992). In the analysis 
of the fluid limit of bandwidth sharing networks a time-scale separation be- 
tween classes usually occurs when one class of traffic reaches equilibrium faster 
than the others, and hence when the fluid limit hits an hyperplane of the state 
space. Simple examples of this phenomenon can be found in (Robert, 2003). A 
more complex example can be found in (Feuillet. 2010). The interesting feature 
in our scaling is the appearance of the averaging principle in the whole state 
space. Similar averaging phenomena have also been studied in statistical physics 
(C. Kipnis, 1991) as well as chemistry and biochemistry (Segel and Slemrod, 
1989) where the kinetics of chemical reactions can be described by systems of 
ordinary differential equations. Usually these works assume that one of the de- 
pendent variable is in steady state with respect to the instantaneous values of the 
other dependent variables. Taking this time-scale separation as an assumption, 
an efficient approximation method called the quasi-steady-state is commonly 
used in that context. This is however in contrast with our situation where we 
show that the time-scale decoupling occurs as a consequence of the scaling of 
the parameters of the transitions of the stochastic processes considered. 

Contribution: 

Our contribution first consists in establishing the convergence in uniformly on 
compact sets for stochastic processes commonly adopted in the modeling and 
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analysis of communications networks under the scaling considered. Since the 
slow part of the processes (surging classes) remains coupled (at a macroscopic 
scale) to the fast part (the remaining classes), such a proof is not standard and 
has to be decomposed in several steps. While preliminary results for monotone 
networks were presented in (Jonckheere, Nunez-Queija and Prabhu, 2010), a 
general proof for general bandwidth sharing networks was still missing. 

Second, we characterize the responses (evolutions of queue length) of differ- 
ent networks to the surge of traffic. We introduce the notion of robust stability, 
which describes a situation when the network can resorb a surge of trafhc by 
eventually reducing the macroscopic state of all classes to 0. We call the set 
of traffic parameters that lead to this condition the robust stability region. We 
characterize this robust stability region for work conserving allocations and for 
monotone allocations. We first show that for work conserving allocations, the 
unstable class, at its macroscopic time scale, "sees" the other classes as hav- 
ing full priority, while the effect of the first class on the other classes gradually 
vanishes (again at a macroscopic time scale) . Hence surging classes tend macro- 
scopically to under the stability condition of the system Pj < 1- The 

situation is more complex for non work-conserving networks, where the behav- 
ior of the unstable class depends in an intricate manner upon that of the other 
classes. In particular, under the usual stability conditions of the network, the 
macroscopic state of surging classes might converge to or to a strictly posi- 
tive number, depending on the conditional distribution of the other classes. For 
monotone networks, we prove that robust stability boils down to stability under 
an allocation giving full priority to stable classes. We illustrate these concepts 
on several simple network topologies. Finally, we use our analytical results to 
build an implementable penalization rule allowing to adapt the level of priority 
of streaming traffic in a network integrating streaming and elastic traffic, such 
as to target a given loss probability threshold/ quality of service. 

The rest of the paper is organized as follows. The model is presented in next 
section. In Section 3, we present the convergence theorem for the considered 
scaling. In Section 4, we analyze the qualitative behavior of networks after a 
traffic surge in different cases and we give numerical examples of applications 
of the main result to bandwidth sharing on some simple network topologies. 
In Section 5, we construct a practical penalization rule for streaming traffic. 
Finally, we conclude in Section 6. 

2. Model 
Notation 

denotes the Zi-norm: 

N 

\x\ = ^ |a-j|. 

i=l 



In the sequel, for x e "L , \ ■ 



imsart-generic ver. 2011/12/01 file: artlcle3.tex date: December 12, 2011 



M. Feuillet et al. /Bandwidth sharing networks with priority scaling 



5 



For x,y G Z , we also use the notation x < y to denote the partial order Xi < yi 
for alH = 1 . . . N. 

2.1. Networks with traffic weighted allocations 

We consider a processor sharing network with N traffic classes. Within each 
of the N traffic classes, resources are shared according to a processor-sharing 
service discipline. The service rates are state-dependent: they may depend on 
the number of flows within the same class, as well as on the numbers of flows 
in all other classes. The service rates of the N traffic classes will be denoted by 
4> = {4>i{'))iLi - Several examples are considered in the next section. Note that the 
service rate function (f) captures the allocation of bandwidth which is determined 
by the specific network topology and congestion control mechanisms. Special 
allocation functions that have received much attention in literature include the 
celebrated max-min fair allocation and the proportional fair allocation. 

We assume that class-i customers arrive subject to a Poisson process of in- 
tensity Xi and require exponentially distributed"'^ service times of mean for 
class-i. The arrival processes of all classes are mutually independent. Our main 
results allow for time- varying arrival rates for the class exhibiting a traffic surge. 
When applicable, we reflect this dependence in the notation by adding the time 
parameter to the arrival rates and then \i{t) is the arrival rate of class-1 at time 
t. For ease of exposition, however, we restrict ourselves to constant arrival rates 
for all classes in this section and will formulate our results with time-varying 
arrival rates in Section 2. 

Let X be the stochastic process describing the number of flows (or calls) in 
progress. In the absence of priority mechanisms, and under the assumptions of 
Poisson arrivals and exponential flow sizes, X is a multi-dimensional birth and 
death process with transition rates: 

q{x,x - e,) = fii4)i{x), 
q{x,x + Ci) = Xi, 

with X £ N^. Assume now that priority mechanisms are employed in the 
network such that the actual bandwidth allocation depends on the variables 
riXi, i ~ \ . . . N rather than simply on Xi, i ~ \ . . . N . Hence, if Xi is thought 
of as a measure of the level of congestion of class-i, a differentiation between 
classes can be enforced by giving different weights to the different classes. (Such 
a differentiation can be enforced at lower time-scales by packet schedulers like 
weighted deficit round robin.) 

It can also be the case that each class of traffic has a limited peak rate 
(because of access constraints for instance). It could then be advantageous for 
providers, in order to meet the demand, to share capacity as a function of the 

-"^Such assumptions are certainly not necessary to obtain the results we are aiming at; 
however, a rigorous generalization would be technically very involved and is beyond the scope 
of the present paper. 
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demanded rates riXi rather than as a function of the number of flows of each 
class in the network. In both configurations, X can now be described as muhi- 
dimensional birth and death process with transition rates: 

qix,x - Ci) = ^ii(f)i{r.x), 
q{x, X + Ci) = Aj, 

where r.x = {riXi)i=i,,,N for some r £ M.^. To avoid confusion, we emphasize 
once more that reflecting the dependence on the control parameters in our 
notation will be more convenient for the purposes in this paper, rather than 
making this dependence implicit, i.e., through the allocation (t>{x) = (j}{r.x). 
The load of class i is given by 

A. 
= — ■ 

p.i 

We shall further suppose that the weights r are chosen for each class propor- 
tionally to the size of the traffic surge, i.e., 

n = ri[\x\) ^ j^. 
Finally, denote the global load of the system by p = X^iLi Pi- 



3. Fluid limits with time scales decoupling 



We model a traffic surge by a large number of initial flows and large arrival rates 
for a subset of classes being (temporarily at least) unstable. Let c < S N. To 
get structural results on the process X, we study the case where: 

1. the number of initial class-i, i = 1 . . .c flows is of order K = \x\. 

2. we scale (accelerate) time by a factor if, 

3. we scale class-z, i = 1 . . .c states by a factor 1/K, 

4. the prioritization weight of class-i, i < c is of order 1/K. 

We now consider a network with several classes of traffic and with class-i, 
i = I . . .c going through a temporary surge of traffic. Recall that we focus on 
a regime where rt = ^ and K ^ oo. We further let denote the (scaled) 
process: 

Y'^it) = (^(^i^),^i...„ (Xf (ift)),=,+i...^.) . (1) 

In the following we show that, as K ^ oo, converges to a stochastic process 
with a deterministic first coordinate, which is a solution of a differential equation 
which we describe in terms of an averaged rate (f>- In the limit, the result implies 
a time-scale separation between the first classes and the other ones. 

Define to be a TV — c dimensional Markov birth-and-death process with 
arrival rates Xi and death rates (j>i{z,-), i = c + 1 . . . N [z G R^) and denote by 
7r^(-) its stationary probability (when it exists). When we do not use a time in- 
dex, we implicitly suppose that we consider stationary versions of the processes. 
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For a given a{t) E R'^, let u{t) G M*"' be the solution (assuming it exists and 
it is unique) of the differential equation: 



Vz = 1 . . . c, Ui{t) 



di{t) - (j),{u{t)), ifui(t)>0, 



0, 



if = 0, 



(2) 



with 0i(z) — X^j^gN"-'^ 4'iiz, y)T^^iy)- To establish our main result, we shall make 
the following assumptions: 

(Al): (f>i{-, Xc+i, ■ • ■ , xn) can be extended to Lipschitz-continuous functions from 
R% \ {0} to M+. 

{A2): for all fixed z, the process is ergodic. We can thus define E*^ the mean 

under the stationary distribution of the process . 
(^3): i/J'*Af(s)ds->a,;(t), i = l...c. 

We can now proceed to state our main result: 

Theorem 3.1. Under the assumptions {Ai), {A2) and {A3), the process ^^(i)^^! 
converges in , uniformly on compact intervals, to the deterministic trajectory 
u{t), i.e.. 



E 



sup \yi^{s) — Ui{s)\ 

0<s<t 



0, K ^ 00, Vi = l. 



(3) 



Moreover, for all times t, and for all bounded continuous functions f: 



lim E sup 

K^oo \ o<s<t 



E^ 



(4) 



The details of the proof are given in the next Section. We underline here the 
main steps: 

• We first prove tightness of the laws of the scaled process, and show that 
the limit-points of are continuous processes. 

• Supposing the convergence in distribution of the first class we characterize 
the limit of the functional ^{{x^{Ks),...,x^{Ks))Gr} ds, VF C and 
prove the limits are unique (and deterministic given the value of the first 
class). A key step is the useful characterization of bimeasures. 

• Finally, we show that converges in distribution towards a deterministic 
process which allows to prove, using the previous step, the convergence in 

, uniformly on compact sets. 
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3.1. Proof of Theorem 3.1 

Step 1: 

For the ease of exposition, we suppose that class- 1 only undergoes a surge of 
traffic. The proof then extends directly to the general case. 

We thus consider the process Y^{t) = j^^^ , {X^ {Kt))i^2...N^ as defined 

by (1). We define N = NU{+cxd} and for each K, we define the following random 
measure on [0,oo) x N^^^: 

/■* 



i/^M(0,t) X r) = / l{(x-(^f,),...,x-(K.))er}ds, VF C and Vi > 0. 

J 

We denote £o(N^~^) the set of measures on [0,00) x N^~^ such that, for all 
measure v in /;o(N^"^) and all i > 0, we have i^((0,t) x W^^) = t. Since N is 
compact, we have that £o(^^^^) is compact and we deduce that {v^ if G N} 
is relatively compact. 

In order to prove the relative compactness of {{Y^ , v^), K e N}, we then 
just have to prove the relative compactness of {Y^ , K £ N}. We define the 
following process 

M^{t) = Y,{t) - - / Ai(s)ds + - / 01 ^Ai,xf (,s) ds (5) 



The martingale characterization of jump processes (see (Rogers and Williams, 
1987)) shows that Mf is a locale martingale and its increasing process is given 

by 

(M,->^ A.Wd. + 5i,f V, (^£^,A-fw)d. 

Using Doob's inequality ^, it follows that il//^ converges in probability to on 
any compact set when K — > 00, i.e., for any T > and any e > 0, 

lim P ( sup \M^{s)\ > £ ) = 0. (6) 



if — >oo 



0<s<t 



We then define Wh the modulus of continuity for any function h defined on [0,t]: 
WhiS) = sup \h{s) - h{u)\. 

s.u<-t; \u — s\<S 



^For any martingale M, using Cauchy Schwartz and Doob's inequality (Darling and Norris, 
2008), we get that; 



E 



sup Ms 

0<a<t 



< E ( sup \Ms\ 

k 0<s<t 



< E ( sup M's 

\0<s<t 

< 4E (Aff ) . 
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Using Equations (5) and (6), we arc able to prove that for any e > and > 0, 
there exists 5 > and A such that for K > A, we have 

The conditions of (Bihingsley, 1999, 7.2 p81) are then fullfilled and the set 
{Y-^ , K 6 N} is relatively compact. Moreover, any limiting point is a continuous 
process. 



Step 2: 

Wc now suppose that (Y^) converges in distribution to a limit Zi. We have 
to characterize any limiting point of the sequence {v^) and then deduce the 
existence and uniqueness of the limit of {v^). In the following, we consider a 
convergent subsequence {Y^\v^'^) and its limit process {Zi,v). 

Let (11, P) be the probability space on which they are defined. We call 
{J^t] the natural filtration of {Z{t), v). We then define 7 such that 

VA eJ;yB<E B{[0, 00)), VC e B{N^-^) -f{AxBxC)^ E{lAiy{B x C)). 

According to (Ethier and Kurtz, 1986, appendix 8), 7 can be extended to a 
measure on J^® S([0,oo)) (g) B(N^~^) and there exists -d such that for all t, 
d{t, .) is a random probability measure on N^~^ and for any B G S(N^~^), 
> 0) is {J'tj-adapted and for any A G J" ® B([0, 00)), 

7(AxB)=E^^ lAisMs,B)dsy (7) 
Meit) = iyi[0,t] X B) - [ ?9(s,B)ds. 



Mb is {7^t}-adapted and continuous. We consider t > s, and D E Fs- We define 
1(7(1^,1*) = l£i(cj)l[s^t) (u) and we have 

¥.{lDv{[s,t) X B)) = 7(1) X [s,t) X B), 
= l{CxB), 



E (^j lc{u)-d{u,B)dv^ , (according to (7)) 
E {^D d{u,B)Au 



Since the previous equality is true for all _D G J^t, it follows that 
E{v{[s,t)xB)\Fs)^E( I d{u,B)du Fs 
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and immediatly, we have 

Then, Mb is a continuous {J-'t}-martingalc. It has finite sample paths and 
then is almost surely identically null. Almost surely, the following equation holds 
for all t, 

VB CN^"\ i/([0,i) X B) = / 'd{s,B)ds. (8) 

Jo 

We have to characterize the random measures 'd{t, .) associated to i/. For any 
uniformly continuous bounded function g on N^^^ and any K G N, we define 

W = {9(^2 m, • ■ • , X'<iKt)) 5(0)) 

[g{X^{Kt),...,Xf{Kt) + e,,...,X^{Kt)) 

z=2 -^0 ^ 

^g{X^{Kt),...,X^{Kt)))As 

-Y.^Ji^ (g {X^iKt), . . . , X^{Kt) - e„ . . . , Xj^{Kt)) 

^g{X^{Kt),...,X§{Kt))) 
0, (r/^(s), (i^t), . . . , X^{Kt)) As. 

As is a martingale, is a martingale. We have that converges 
in distribution to 0. \Ki\-^{g{X^' {Kit), . . . , X^' {Kit)) - g(0)) also converges 
to because g is bounded. As a consequence, the following term 

V A, /* (g (Xf ' {Kt), . . . , ' {Kit) + e„ . . . , X^' {Kit)) 
i=2 -'o V ^ ^ 

-g(x^'{Kit),...,X'^^{Kit))^ds 

-g(x^'{Kit),...,X^^{Kit)) 
0, (Fi^' (s), ' {Kit), ...,X^^ {Kit)) ds 
also converges in distribution to 0. But, by the continuous mapping theorem 
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and (8), it converges in distribution to 



t N / 

H r"'' 9{y + et)-9{y) 



i=2 \ ygN"-i 



+ Y (5(2/ - eO - 5(2/))0« (^i('S),y) n?(s,2/)ds. 



Consequently, this is null almost surely for all t and wc have then, for Lebesgue- 
almost every i. 



where is the infinitesimal generator of {U^^^^\s)). This proves exactly 

that .) is invariant for f/^^'*'. By uniqueness of the invariant distribution 
of ({/^^^*-'(s)), this implies that, given Zi, §{1, ■) is a deterministic measure for 
all t. We can deduce that, if iY^'^) is a converging subsequence, then (w^') is 
also converging and its limit is a random measure in £o(N^~^). This implies in 
particular that (w^') is tight in £o(N^^^). We can now proceed of the last part 
of this step. 

We consider £ > 0, 77 > and t>{). Because the sequence {v^^) is tight in 
£o(r^^^^)i there exists k > and a compact F C N^^^ such that: 



Because Zi is almost surely continuous and / is Lipschitz-continuous, we 
have 




We deduce immediately that 





\ i>K,t^er,s>t 
Since / is bounded, we can deduce: 



P sup >e <7;/2. 






f{Z^{s),y)u{dsxdv) >2e\\f\\ < 77. 
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According to (8), there exists a family {"dit, .)) of random measures on ^ 
such that 



sup 

0<s<t 



f f{Y^^{u),Xf^{Kiu))- nZi{u),y)d{u,y)du 



converges in probabihty to when Ki tends to infinity. 

Since / is bounded, we can apply the dominated convergence theorem and 
we have that 



K, 







sup 




\ 0<s<i 





We further have that 



lim E sup 

X->oo \o<s<t 



Step 2 is complete. 



/(n^'(u),Xf'(A») 
^ f {Zi{u),y)'d{u,y)du 



0. 



du 



= 0. 



Step 3: 

Using the martingale decomposition of 
Xi<iKt) 



1 r^* 

xi+MK{t) + - / Af(s)d, 
^ Jq 

X^{s) 



K 



1 



,...,Xf(,s) d,s 



As already remarked in step 1, since (j) is bounded it follows that 

which implies using Doob's inequality that there exists a constant A' such that 
for K big enough: 



E 



sup Mx{s) 

0<s<t 



< A'\ — < e. 

- y K - 



Using the convergence of the arrival process together with the convergence of 
the martingale Mk, we obtain the uniform integrability of . (The tightness of 
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has already been obtained in step 1). Now consider a converging subsequence 
Y^' towards Zi. Using the results of step 2, the convergence of the arrival 
process together and the convergence of the martingale Mk, we obtain that Zi 
must satisfy: 

Zi{t)^xi+0 + ai{t)- [ 0i(Zi(s))ds. 



Hence the limit is unique and deterministic. 

This in turn shows the convergence of Y^ in distribution and completely 
characterize the measure •& introduced in Step 2 as a deterministic measure. We 
can now prove the convergence in L^. Let e be given. Define the error estimate: 

nxit) ^ sup \Y^^{s) - ui{s)\. 

0<s<t 



Define now the noise amplitude as: 

MKit) = sup \MKis)\ . 
a<s<t 

Using the convergence of the intensity of the arrival process, 
riKit) <MK{t)+e 



sup 

s<t 



1 



K 



dz 



01 (u(z)) dz 



Using step 2, (pi being Lipschitz and bounded, for K large enough: 



E sup 



Ks 



Xli_z) 
K 



,Xf{z) ]dz- I (j)i{u{z))dz 



which concludes the proof for the convergence of Y]^ 



K 



4. Qualitative behaviors of the limiting processes 

We describe thereafter the different qualitative behaviors that may occur de- 
pending on the traffic conditions. Assume that class-i, i = 1 . . . , c have en- 
tered a traffic surge . Under the scaling considered in Theorem 3.1, we observe 
three qualitative types of behaviors for the network responses, which are com- 
pletely characterized using the stationary distributions of the family of processes 
U^,i = 2,...,N. Defining 

Vx e M^", Vi G {1, . . .,K}, S^{x) = A, - ^l^Mx), (9) 

let A the set of positive solutions of the equation 

S{x) ^ 0. 
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Given the classical results on asymptotic stability of non-linear autonomous 
systems, we can partially classify the possible situations using in particular the 
Hartman-Grobman theorem (see for instance (Hartman, I960)). For a flow 
6, we write DS{x) < if the linearization of 6 has only eigenvalues with strictly 
negative real parts and no eigenvalue on the unit complex circle. Assume in the 
following 5 is C^. We have the possible behaviors: 

1. The network continues to see class-i, i > c saturated (at a macroscopic 
time and space scales) , even after any large (macroscopic) amount of time. 
A sufficient condition is that there exists x £ A such that D5{x) < and 
a; > 0, with initial conditions sufficiently close to x. Then the differential 
equation is asymptotically stable with stable point x > 0. In this case, 
limits in time and K cannot commute since there is always a part of the 
bandwidth of the network used for surging classes, while taking first the 
limit in time and then the limit in K always converge to the system with 
allocation 0i(O, . . ■ , 0, Xc+i, ■ ■ ■ , xn) for stable classes. 

2. The differential equation (2) governing the dynamic of u is unstable, which 
means that the traffic surge cannot be absorbed and keeps building up. It 
might lead to the instability of other classes in the network. 

3. The traffic surge will be absorbed at macroscopic time, i.e., the differential 
equation is asymptotically stable with stable point 0. Necessary conditions 
for this situation are that: 

(a) OeA 

(b) D6{0) < 0. 

(c) the initial condition is close enough to 0, or ^ = {0}. 

In this case, note that the stationary measure of {Y/^ {t))iyc converges 
when ii' oo to the stationary measure of the original system with al- 
location (^i(0, . . . , 0, Xc+i, . . . , Xn)), which boils down to the fact that the 
limit in time and in K commute for classes c 4- 1 to N . 

4. The traffic surge will be absorbed at macroscopic time, i.e., there exists a 
T > 0, such that for all t > T, the solution is of the differential equation 
is null. 

4.I. Robust bandwidth sharing networks 

Bandwidth sharing networks constitute a natural extension of a multi-class pro- 
cessor sharing queue, and have become a standard stochastic model for the 
flow level dynamics of Internet congestion control (they were introduced by 
(Massouhe and Roberts, 2002)). 

Consider for example the tree network represented on the left of Figure 1, 
with two traffic routes, each passing through a dedicated link, followed by a 
common link. If each dedicated link has a capacity Ci < 1, i — I, 2, and the 
common link has capacity 1, the flow on each route gets a capacity (j)i{x) that 
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lies in the polyhedron C: 



1 



(10) 



0j(a;) < Ci, i = 1,2. 



(11) 



Another example of interest is the linear network represented on the right of 
Figure 1 with 3 routes sharing two links. While the first route passes through 
both links, routes 2 and 3 only use one of the links (one eaeh). This gives the 
following capacity constraints: 



In general, like for the specific foregoing examples, the capacity constraints 
determine the space over which a network controller can choose a desired allo- 
cation function. It has been argued by (Kelly, MauUoo and Tan. 1998) that a 
good approximation of current congestion control algorithms such as TCP (the 
Internet's predominant protocol for controlling congestion) can be obtained by 
using the weighted proportional fair allocation, which solves an optimization 
problem for each vector x of instantaneous numbers of flows. Specifically, the 
weighted proportional fair allocation 77 (x) for state vector x maximizes 



where the weights Wi are class-dependent control parameters. 

Remark 4.1. By definition of this optimization program, if 4>{') ~ vi') 
standard (unweighed) proportional fair allocation with Wi = 1, then the alloca- 
tion (jf{x) = 4>{r.x) corresponds to the weighted proportional fair allocation with 
weights Wi = ri. 

This framework has been generalized to so-called weighted a-fair allocations, 
which provide flexibility to model different levels of fairness in the network. 
Another important alternative is the balanced fair allocation (Bonald et al., 
2006), which allows a closed form expression for the stationary distribution 
of the numbers of flows in progress. In addition, the balanced fair allocation 
gives a good approximation of the proportional fair allocation while being easily 
evaluated, which is attractive for performance evaluation. 



'^1(2:) + <t>2{x) < Ci, 

^i(a;) + (/)3(a;) < C2. 



(12) 
(13) 



N 




Fig 1. Tree network and linear network 
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Remind that all a-fair bandwidth sharing arc stable for a > (in the 
sense that the process X is positive recurrent) (Bonald and Massoulie, 2001; 
de Veciana, Konstantopoulos and Lee, 2001) if 

peS = {v,A7^<C}. 

We now refine the concept of stability by saying that the network is sta- 
ble if classes undergoing a surge eventually drain while the other classes stay 
stochastically stable. More formally, we say that 

Definition 4.1. The network is robust stable if: 

E-[X-{\x\t),] 
iim sup lim sup j — j = U. 

t-i-oo |x|— s-oo \x\ 

Define further the robust stability region as the set of parameters such that 
the network is robust stable, i.e.: 

5'- = L G : hmsup sup ^ q] . (14) 

In the sequel, we first show that for work-conserving allocations, the robust 
stability set coincides with the usual stability set (except possibly on a negligible 
set of parameters). 

On the other hand, the situation is much more complex for non-work-conserving 
allocations where even with appropriate weighted allocations, surges might not 
be asymptotically transparent to the other classes while classes undergoing 
surges might get asymptotically strictly more bandwidth than in the case of 
the allocation giving full priority to stable classes. Wc give an example of this 
phenomenon later on. 

For monotonic networks, wc however prove that the robust stability coincide 
with the set of parameters under which a "priority allocation" is stable. 

4-2. Existence of the fluid limit 

Lemma 4.1. If p E S, then the processes are positive recurrent for any z, 
and Theorem 3.1 applies. 

Proof. Asymptotically, when — > oo, i > c the allocation allocated to classes 
i > c coincides with the allocation 0(0, . . . , 0, Xc+i, . . . , xn) which is stable under 
the usual conditions of traffic. □ 

4-3. Work-conserving allocations 

Consider a work conserving allocation such that 

N 

VXT^O, ^4>^{x) = l. 
1=1 
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Every work-conserving allocation has the same stability region, namely 

N 

'i=i 

If the (usual) stability condition is satisfied, then the priority mechanism con- 
sidered is asymptotically equivalent to giving full priority to class-i, i > c-l- 1. 
In other words, the fluid limit obtained for class-1 is in that case the same as 
the fluid limit of an allocation that gives a full priority to class-i, i > c 4- 1, 
which we prove in the following Proposition. 

Proposition 4.1. For a work conserving network, S'' — S (except possibly on 
the frontier of the stability sets) and: 

y,^(t) u,{t) = |^u,(o) + J2 ^ 

Proof. Assume X^i^i Pi < 1 in which case the network is stable. Fix zi e M. 
Using the conservation of the rates at equilibrium for the process U^^ ( which 
boils down in the Markovian context to saying that at equilibrium the drift of 
y — > 2/i should be 0), we can write that: 

N N 
2— c+l y i—c 

We now calculate 4>i for zi > 0: 

(l^iix) =^(f>i {zi , 2/)7r'^i (y) , 
y 

,^,(x)=^(l- ^ ^,(zi,y))7r^Hy), 
y j>c+i 

j>c+l 

Hence, the capacity seen asymptotically by class-1 is (1 — J2j>c+i Pj)^ which 
concludes the proof. 




Example: one link v^rith the DPS allocation The simplest instance of a 
network consists of one link shared by several competing classes of traffic. If the 
initial policy is supposed to be the classical processor sharing policy: 4>i{x) = 
then the prioritized version of the model becomes the so-called discriminatory 
processor sharing (DPS): (t)i{r.x) = ^'''^'^ . 

Consider a single link of capacity 1 shared by three classes. The bandwidth is 
allocated according to DPS with weight for class i, i = 1, 2, 3. Proposition 4.1 
says that ui{t) is a straight line with slope Ai — ^i(l — {p2 + Ps))- This behavior 
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Fig 2. DPS with three classes: scaling of class-1 (left) and of class-2 (right) 




is illustrated in Figure 2, for which Ai = 0.5, /ii = 1,P2 = 0.3, /33 = 0.1. The 
slope calculated using the proposition is thus 0.1, which is verified in the figure. 

In Figure 2, we plot the empirical mean of class-2 at a macroscopic scale, (i.e. 
i /t fO^^W) d^) for a temporal window of s = 0.1. 

4.4- Monotone allocations 

Define the allocation ip giving full priority to class c + 1 to TV, given by 

ip{x) = (^(0, . . . , 0, Xc+i, . . . , xat), if > 0, for some i > c. 
ipix) = 4>{x), otherwise. 

Denote S{il]) the stability region of the network with allocation t/). 

Proposition 4.2. Consider a monotonic allocation (i.e. such that is decreas- 
ing in Xj, j i). 

If (5(0) < where S{x) is defined by (9), then the network is robust stable 
and surging classes do not influence asymptotically stable classes. Conversely if 
6{0) > 0, the network is not robust stable. 

Moreover S{0) < if the network with allocation tp is stable i.e.: 

Proof. Using stochastic comparisons (see (Borst, Jonckheere and Leskela, 2008) 
for more details on stochastic comparisons of multidimensional birth-and-dcath 
processes with monotonic allocations), we obtain that 

f/° <.t c/^vz, 

which implies that Vz, 4>i{z) > 4>i{0), This in turn implies that 

^Mt) < ^(0) < 0, Vt > 0, such that u,{t) > 0. 

This implies that Ui will reach in finite time. 

The reverse statement follows along the same lines. □ 



imsart-generic ver. 2011/12/01 file: article3.tex date: December 12, 2011 



M. Feuillet et al. /Bandwidth sharing networks with priority scaling 



19 



Example: a tree network Let us consider the tree network shown in Figure 1 
with ci = 0.4 and C2 = 0.8. We shall assume the following bandwidth allocation: 
Define Si = {(xi,X2) : {rixi + r2X2)ci < rixi}. For a;i > and X2 > 0, 

and (/)2 = 1 — 01. 

For this network, the allocation becomes a strict priority allocation for class- 
2 when ri = 0, in which case class- 1 gets capacity ci if there arc no class-2 
flows, and 1 — C2 otherwise. Thus, for a fixed value of p2, class-1 is stable if 
pi < 1^1 — |^^ci + ^(l — C2). The stability regions for ri = and ri > are 
shown in Figure 3. 



Fig 3. Partitioning of the stability region for the tree network 



P2 




The dynamics of ui{t) for two different values of pi - one in each region - is 
plotted in Figure 4, for which p2 = 0.5. 



Fig 4. Tree network: scaling of class-l (left) and of class-2 (right) 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 



For class 2, when the priority allocation is stable the dynamics of the average 
number of customers converges to the one of the priority allocation, that is 
P2/(c2 — P2), as is illustrated in Figure 4. 
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In Figure 5, wc show how class- 1 is actually favored by asymptotically using 
the bandwidth of class-2, compared to the case where class-2 is given a strict 
priority. 

Fig 5. Tree network: comparisons of trajectories of class-l for a proportional fair allocation 
and a priority (to class-2) allocation 




' ' 1 ' 1 ' 1 1 1 

5 10 15 20 25 30 35 40 45 
t 



4-5. Non-monotone networks 

For non-monotonic networks, some unusual behaviors can be observed; for in- 
stance the usual fluid limit may exhibit the following behavior: one class of traffic 
can reach at the fluid scale and stay at for a finite time before increasing 
again. We here show on an example how the priority scaling avoids this kind of 
behavior. 

Here, we consider a linear network with two links and three classes of flows 
as shown in Figure 1. 

Let ai be the capacity allocated to a flow of class-i. The capacity allocated to 
class-i, 4>i, is then Xtai. The bandwidth is allocated according to the weighted 
proportional fair allocation, that is, {ai,a2,a:^) is the solution of the following 
maximization problem: 

maximize rixi log(ai) -I- r2X2 log(a2) -I- 7-3X3 log(Q;3) 
subject to xiai + X2ce2 < ci, 
xiai + xsas < C2, 

where is the weight of class i. 

Consider the following network parameters with arrival rates and service 
rates: 

ci = C2 = 1, ri = r2 = rs = 1, 
Xi(0) = 10-K,X2{0)=K,X3{0)^K, 
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that is, the three classes are unstable at the beginning. The trajectories of the 
number of flows as a function of the scaled time for K = 10000 is shown in 
Figure 6. 



Fig 6. Usual fluid (left) and priority scaling (right) for the linear network. 




Concerning the usual fluid limit, we observe that Class 2 becomes stable 
around the 2 time unit mark. It then becomes unstable around the 15 time unit 
mark, and becomes stable again around the 40 time unit mark. This behavior 
can be explained as follows. 

In the first part of the trajectory, class 2 gets sufficient capacity to drain out 
while the number of class 3 flows grows. In the second part, since Link 2 is not 
a bottleneck for Class 3 flows and Class 1 gets smaller. Class 3 gets a larger 
share on Link 2. Meanwhile, the arrival and service rates of Class 2 and Class 
3 being the same, the imbalance in the rate allocation means that Class 2 is 
now unstable and its number of flows starts to grow until it reaches the same 
number as that of Class 3 flows at which time they both share the link capacity 
equally. Since the network is stable, all the three classes drain out eventually. 

Due to the priority mechanism employed to penalize unstable classes, such a 
phenomenon does not happen in the case of the priority scaling as illustrated 
in the right-hand Figure. 

5. Integration of streaming and elastic traffic 

Consider now a system where two intrinsically different types of traffic - "stream- 
ing" and "elastic" traffic - coexist and share a given link. Such models have been 
considered by (Nunez-Queija, van den Berg and Mandjes, 1999; Delcoigne, Proutiere and Regnie, 
2004; Bonald and Proutiere, 2004). It is natural to equip streaming traffic with 
a fixed required rate, say, c per flow. Giving priority to streaming traffic (class 
2) the allocation of service may be chosen as: 

01 (.t) = max I ^— i , 1 — CX2 ] , 

\rixi + CX2 ) 

4'2{x) = CX2, 
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where the parameter ri quantifies the level of priority. The alloeated eapacity 
cannot exceed the total capacity. If the latter is normalized to 1, the state space 
must be restricted to states X2 such that 

Mx) + Mx) < 1- 

Then, if the number of current streaming flows X2 is such that 4'i{x + e-i) + 
(^2(3^ + 62) > 1, arriving streaming flows must be blocked from the network. 

In this example, if accepted in the network, the capacity allocation of class- 
2 flows is CX2 independently of the number of flows of class- 1, for all values of 
ri > 0, whereas the capacity allocation of class- 1 flows depends on the number of 
class-2 flows such that 01 (xi, x^) = rxXxjirxXx-Vcx^)- However, class-2 flows are 
admitted only if there is sufficient capacity, that is, if ;;jY(t]T^lfJ7+T) 

(a;2 + l) < 1. 

Denote 

Szx ^\x2^■ T h CX2 < 1 

1^ Zi + CX2 

the state space of class-2 conditioned on u\{i) ~ z\. Define p2 = The 
process J7|^ is birth-death process with birth rate A2 and death rate [i^cx^-, and 
whose stationary distribution is given by 

For the priority allocation, class-1 is stable if and only if p\ < 7r2(0). Thus, if 
Pi < ''■2(0), then the limit point of ui{t) is 0, and if 7r2(0) < pi < 1, then the 
limit point is positive. 

Performing the scaling previously defined, remark that the state space de- 
pends for a fixed macroscopic state zi on both zi and c. We can apply Theorem 
3.1 with (j>i being defined by: 

Mz)- E -^^^c{z,i 

^ Zi+CX2X2\ 

where C{zi) = {J2x2£S^ ^^^^ ^^^^ ^^^^ small (c << 1), we 

might consider as a reasonable approximation a Poisson distribution for class-2, 
whatever the state of class-1. In that case, ^1 takes a slightly simpler form. After 
simple calculations: 

^:i/tf'M^^~^exp(M) du 

czi) = H{zi) = y-^— — . 

P2 exp(p2) 

This allows a recursive evaluation for integer- valued zi. Using simple calculus, 
for n e N: 

77 + 1 

H{n + 1) = ^—{l-Hin)) 

P2 
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We can also evaluate H in terms of special functions: 

n\ — nT{n, —1) 



Hin) = 



(-P2)"exp(p2)' 



where T{n, —1) is the incomplete F function. 

In Figure 7, we plot the ui{t) for pi — 0.6, p2 — 0.2. and c = 0.01. 



Fig 7. Streaming and elastic traffic: scaling of the elastic traffic 
1 

Pl=0.6,Jf = 1000 
P1=0.6,K = 10000 
"■9 n pi=0.6,K = 




10 12 14 



5. J. Quality- of- Service guarantee 

The Quality-of-Service for streaming flows is mainly characterized by the prob- 
ability that an incoming flow does not find sufficient capacity in the network 
for the fiow to be accepted. In networks in which streaming and elastic traffic 
do not interact, this probability can be computed using an Erlang Fixed-Point 
approximation (Kelly, 1986). However, in the context of the present example, 
the interaction of these two types of traffic makes it more difficult to apply these 
fixed-point approximations, mainly due to the fact that the state space of the 
elastic fiows is unbounded. However, in the limiting regime under consideration, 
we can come up with a rule-of-thumb that can be used to guarantee a blocking 
probability smaller than a desired value. 

First, we consider a single link whose capacity is shared by the two types of 
flows. Let pm denote the desired maximal blocking probability of class-2 flows. 
We shall set the priority level of class- 1 (by varying ri) such that the probability 
of blocking of class-2 is always less than p„i . 

For ui(t) = zi, an arrival of class-2 is blocked if and only if — -I- 

c{x2 + 1) < 1, which is equivalent to - 1 < ^2 < ^-7^. The term is 
the number of circuits of size c available when the total capacity is 1 — zi. Thus, 
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is the maximum number of simultaneous flows of class-2 in the system, and an 
arrival of class-2 is blocked if and only if the number of flows of class-2 is N2- 

Let g{n) denote the blocking probability when the number of circuits is the 
network is n. From the Erlang-B formula, 

g{n^ - 



The inverse function g^^ipm) gives the minimum number of circuits required to 
ensure a blocking probability smaller than pm- In order to guarantee a maximal 
blocking of the number of circuits, N2 has to be larger than g~^{pm) at all 
instant of time, which leads us to the following necessary and sufficient condition 
for guaranteeing the Quality-of-Scrvicc of class-2 flows: 

ui := sup ui{t) < 1 - c\g^^{p„i)]. 

0<t<oo 

We can ensure the above inequality by scaling the process ui{t) by a factor 
"'""'^^^-^ (P"*)] ^ This, in turn, can be achieved by scaling the priority level (or, 
equivalently, ri) by this very same factor. This additional scaling results in 
a larger share of the bandwidth for class-1 flows in case 1 — c\g~^{pm)'] > 
Ml. Conversly, if 1 — c\g^^{p,n)^ < "i, the priority level of class-1 flows is 
appropriately decreased so that the blocking probability constraint of class-2 
flows is not violated. 

Using the monotonicity of 0i in its first variable, we get that if Ai > (ui (0)), 
then ui{t) converges monotonically to its limit point. Hence, 



Ml = 




if Al <<^i(ui(0)); 
otherwise. 



Remark 5.1. The blocking probability for a given value of zi is in fact a con- 
ditional blocking probability in the sense that it is the fraction of calls dropped 
when the class-1 flows take away a capacity of zi. The unconditional blocking 
probability of class-2 flows can be computed by integrating over zi , which is 
rather conservative. An alternative scaling could be constructed such that only 
the unconditional blocking probability satisfies a given constraint. 

Remark 5.2. In a network of links shared by several classes of streaming flows 
and one class of elastic flow, we could use fixed-point approximations to compute 
the blocking probability for the different classes of streaming flows as a function 
of zi. Assuming that this probability is increasing in zi, we could then compute 
the maximum value that zi can attain without the streaming classes violating 
their individual blocking probability. 

6. Conclusions 

We analyzed the flow-level performance of multi-class communication networks 
when one of the classes undergoes a traffic surge. We showed that, under an 
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appropriate scaling of space and time, the dynamics of the temporarily unstable 
class can be described by a deterministic differential equation in which the time 
derivative at a given point depends on the conditional stationary distribution of 
the other classes calculated at that point. For work-conserving allocations, the 
differential equation is the same as the one of the network in which other classes 
have strict priority over the temporarily unstable class, that is, the scaled process 
evolves linearly and is either absorbed at zero or grows indefinitely depending 
on whether the network is stable or not. 

For non-work conserving allocations, the trajectory is much more complex to 
describe as it depends on the mean residual bandwidth left over by the other 
classes which in turn depends on the current state of the first class. The limit 
point of the fluid trajectory can hence be non-zero and finite. We characterized 
the robust stability region of monotone allocations. We illustrated this behavior 
through several examples of network topologies and bandwidth allocations that 
are commonly used to model communication networks. 

The time-space-transitions scaling that we considered raises several open 
questions which would give a better understanding of the network dynamics. 
In particular, finding necessary and sufficient conditions for the limit point of 
non work-conserving allocations to be zero would constitute a very interesting 
result. Also, error bounds estimates would be necessary to obtain a reliable 
performance evaluation tool. 
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